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We reanalyze a new quintessence scenario in a brane world model, assuming that a quintessence 
scalar field is confined in our 3-dimensional brane world. We study three typical quintessence models 
: (1) an inverse-power- law potential, (2) an exponential potential, and (3) kinetic-term quintessence 
(fc-essence) model. With an inverse power law potential model (V((j>) = fi a+4 cj)~ a ), we show that in 
the quadratic dominant stage, the density parameter of a scalar field 0,$ decreases as a -4 ( a_2 )/( a + 2 ) 
for 2 < a < 6, which is followed by the conventional quintessence scenario. This feature provides 
us wider initial conditions for a successful quintessence. In fact, even if the universe is initially in a 
scalar-field dominant, it eventually evolves into a radiation dominant era in the p 2 -dominant stage. 
Assuming an equipartition condition, we discuss constraints on parameters, resulting that a > 4 
is required. This constraint also restricts the value of the 5-dimensional Planck mass, e.g. 4 x 
10 _14 m4<m.5<3 x 10~ 13 m4 for a = 5. For an exponential potential model V = fj 4 exp(— X<j}/rrn), 
we may not find a natural and successful quintessence scenario as it is. While, for a kinetic-term 
quintessence, we find a tracking solution even in /9 2 -dominant stage, rather than the fi^-decreasing 
solution for an inverse-power-law potential. Then we do find a little advantage in a brane world. 
Only the density parameter increases more slowly in the p 2 -dominant stage, which provides a wider 
initial condition for a successful quintessence. 
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I. INTRODUCTION 

Recent observation of the angular spectrum of the cos- 
mic microwave background (CMB) in a wide range of t, 
besides the results for small £ obtained by COBE, sug- 
gests a flat universe By combining this result with 
the measurements of high-redshift supernovae [g] , which 
predicts that the expansion of the universe is accelerating 
now, we are forced to recognize the existence of a cosmo- 
logical constant, or a kind of dark energy, which value is 
almost the same order of magnitude as the present mass 
density of the universe. From the viewpoint of particle 
physics, however, it is quite difficult to explain such a tiny 
value by 14 orders of magnitude smaller than the elec- 
troweak scale. The failure of theoretical explanations for 
the present value of the cosmological constant is known 
as the so-called "Cosmological Constant Problem." j| 

Since there is no plausible theoretical idea to explain 
such a small cosmological constant, it may be more plau- 
sible that such a tiny value is achieved through the dy- 
namics of a fundamental field. This idea is called a decay- 
ing cosmological constant 0, ||. Among such models, 
the so-called quintessence seems to be more promising. It 
shows an interesting property, i.e. a quintessence field fol- 
lows a common evolutionary track just as an attractor in 
a wide range of the initial conditions, so that the cosmol- 
ogy is extremely insensitive to the initial conditions 
§,§|,§. Some of them (so-called "tracker-fields") may 
avoid the coincidence problem. According to a successful 



quintessence scenario, the energy of a scalar field tracks 
the radiation energy (or matter energy) for a rather long 
time in order not to affect the nucleosynthesis at the ra- 
diation dominant era, and the structure formation at the 
matter dominant era, and then becomes dominant just 
before the present time. 

There seems to be, however, still a kind of fine-tuning 
problem for these models. Although the initial energy 
of a scalar field could be the same as that of radiation 
fluid, most contribution is from its kinetic energy. In the 
kinetic dominant case, the scalar field behaves as a mass- 
less field, which is equivalent to a stiff matter. Hence its 
energy density drops as a -6 which is much faster than 
the radiation energy. Then the kinetic term drops in 
the radiation dominant universe until a tracking solu- 
tion is found. This is why the density parameter of a 
scalar field decreases before reaching a tracking solution. 
Since we have to tune the mass scale of a potential for a 
successful quintessence, the potential term cannot be so 
large. Hence, an equipartition condition, which may be 
expected in the early stage of the universe, seems to be 
unlikely in such a model. Another unsatisfactory point is 
that if the universe starts in a scalar-field dominant con- 
dition, the radiation dominant universe is never recov- 
ered. Then we are not able to find the present universe. 
Some modified models within the conventional gravity 
theory have been proposed to solve such problems p0[ , 

& 

As for the early stage of the universe, recently we have 
a new interesting idea, which is called brane cosmology. 
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In a brane world scenario, our universe is embedded in 
higher dimensions and standard-model particles are con- 
fined to four-dimensional hypersurfaces (3-branes), while 
gravity is propagating in higher-dimensions (a bulk) fl2| | 
- JTq] . Among them Randall-Sundrum's second model is 
very interesting because it provides a new type of com- 
pactification of gravity. Assuming the 3-brane has pos- 
itive tension and is embedded in 5-dimensional anti de- 
Sitter bulk spacetime, the conventional four-dimensional 
gravity theory is recovered, even though the extra di- 
mension is not compact |Q. While, in a high energy 
region, gravity theory is very different from the four di- 



mensional Einstein theory 
the geom etrical aspects anc 



16[| . Many authors discussed 
its dynamics (For a review, 



> dynamics (rt 

see, |7|), as well as cosmology |J, II, & 

The main modified point from conventional cosmol- 
ogy is the appearance of the quadratic term of energy- 
momentum and dark radiation. Since the quadratic term 
changes the dynamics of the universe in its early stage, 
we expect some improvement for a quintessence scenario. 
In pi], it is shown that the quadratic term indeed dras- 
tically changes the evolution of a scalar field with an 
inverse-power-law potential and then the density param- 
eter of a scalar field decreases in time until conventional 
cosmology is recovered. This result helps to construct a 
more natural and successful quintessence scenario. 

The purpose of the present paper is to present a full 
analysis of the previous letter, including a numerical 
analysis, and to study other types of quintessence mod- 
els in the brane world scenario. In Secjn], we present 
the basic equations assuming the Randall-Sundrum II 
model. We then study three typical quintessence mod- 
els, an inverse-power-law potential model (in Sec |IIip , an 
exponential potential model (in Sec. IV), and a kinetic 
term quintessence (fc-essence) model (in Sec.[y|). Sec. VI 
is devoted to the conclusions. 



II. BASIC EQUATIONS 

We start with the Randall-Sundrum type II brane sce- 
nario because the model is simple and concrete. It 
is, however, worthwhile noting that the present mecha- 
nism may also work in other types of brane world mod- 
els, in which a quadratic term of energy-momentum ten- 
sor generically appears. In the brane world, all matter 
fields and forces except gravity are confined on the 3- 
brane in a higher-dimensional spacetime. As for gravity, 
in contrast to the other type of brane scenario, in the 
Randall-Sundrum type II model, the extra-dimension is 
not compactified, but gravity is confined in the brane, 
showing the Newtonian gravity in our world. Since the 
gravity is confined in the brane, it can be described by 
the intrinsic metric of the brane spacetime. By use of 
Israel's thin shell formalism and assuming Zi symmetry, 
the gravitational equations on the 3-brane is given by 

(4) G>„ = - (4) A 5ai „ + k\T^ + - E^, (2.1) 



where ^G^ v is the Einstein tensor with respect to the in- 
trinsic metric g^, is the 4-dimensional cosmological 
constant, T Mt , represents the energy- momentum tensor of 
matter fields confined on the brane and 7r M „ is quadratic 
in T^ v |l(J. is a part of the 5-dimensional Weyl 

tensor and carries some information about bulk geom- 
etry. k 4 2 = 87rG 4 and k 5 2 = are 4-dimensional 
and 5-dimensional gravitational constants, respectively. 
In what follows, we use the 4-dimensional Planck mass 

TO4 = k 4 _1 = (2.4 x 10 18 GeV) and the 5-dimensional 

2/3 

Planck mass 11J5 = k 5 , which could be much smaller 
than 7714. We also assume that ^ 4 ^A vanishes. 

Assuming the Friedmann-Robertson- Walker spacetime 
in o ur b rane world, we find the Friedmann equations from 
Eq.(O) as 
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where a is a scale factor of the Universe, H = a/ a is its 
Hubble parameter, k is a curvature constant, P and p 
are the total pressure and total energy density of matter 
fields. C is a constant, which describes "dark" radiation 
coming from E^ v [|l8|. In what follows, we consider only 
the flat Friedmann model (k = 0) for simplicity. 

As for matter fields on the brane, we consider a scalar 
field 4> as well as the conventional radiation and matter 
fluids, i.e. p — p^ + p T + p m , where p^, p T and p m are the 
energy densities of scalar field 0, of radiation fluid and 
of matter one, respectively. Although we can consider 
a 5-dimensional scalar field living in the bulk |2^| , we 
shall focus only on a 4-dimensional scalar field confined 
on the brane. The origin of such a scalar field might be 
a condensation of some fermions confined on the brane. 

Since the energy of each field on the brane is conserved 
in the present model, we find the dynamical equation for 
such a 4-dimensional scalar field as a conventional one, 



• dV 

ZHcj) + — = 0, 



(2.4) 



where V is a potential of the scalar field. The energy 
density of the scalar field is 

Pt = l<P 2 + V(4>), (2.5) 

and for the radiation and matter fluids, we have 

p\ + AH Px = 0, (2.6) 
Pm + 3H Pm = 0. (2.7) 

As the Universe expands, the energy density decreases. 
This means that the quadratic term was very important 
in the early stage of the Universe. Comparing two terms 
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(the conventional energy density term and the quadratic 
one), we find that the quadratic term dominates when 



p > p c — 12m 5 6 /m 2 . 



(2.8) 



When the quadratic term is dominant, the expansion law 
of the Universe is modified. For example, in the radiation 
dominant era, i.e. p t ^> p m ,p$, the Universe expands 
as a a i 1 / 4 , in contrast with i 1 / 2 in the conventional 
radiation dominant era. 

Since we are interested in a quintessence scenario and 
then the dynamical behavior of a scalar field, we shall 
calculate the density parameter of the scalar field, which 
is a ratio of the energy density of the scalar field to that 
of total energy density (fi^,). Ignoring dark radiation, we 
find 



Hp, 
P 2 



(2.9) 



When the quadratic term dominates in the early stage, 
the Hubble expansion rate decreases. As a result, the 
friction term in Eq. (2.4) becomes small and then the dy- 
namics of the scalar field will drastically change from the 
conventional one even for the same potential. It is easily 
expected that a kinetic term will play a more important 
role in the energy density of a scalar field. Consequently, 
if (f) 2 > 4V(0), the density parameter of a scalar field 
will decrease in time. In particular, this feature may be 
very important in a quintessence scenario, because most 
quintessence models assume very small potential energy, 
initially, for a successful scenario. We will show that 
this interesting feature of a scalar field is found in some 
quintessence models and the initial conditions for a suc- 
cessful scenario become much wider. 

In the following sections, assuming that a quintessence 
field exists in a quadratic term dominant era, we inves- 
tigate the dynamical behavior of a quintessence scalar 
field and search for a natural and successful quintessence 
scenario. We study three types of quintessence fields: 
models with an inverse power law potential, with an ex- 
ponential potential and a kinetically driven model. 



III. INVERSE POWER LAW POTENTIAL 

In this section, we investigate a model with an inverse 
power law potential @, i.e. 



(3.1) 



where p is a typical mass scale of the potential, which is 
not fixed here. Although this potential is not renormaliz- 
able, it may appear as an effective potential for a kind of 
fermion condensation in a supersymmetric QCD model 
p3[ . The SU(N C ) gauge symmetry is broken by a pair 
condensation of A^j-flavor quarks, and the effective po- 
tential for a fermion condensate field 4> is given by (3.1) 
with a = 2{N C + N f )/(N C - N f ). 



The quintessence scenario with this potential has been 
well studied in the conventional universe ||. We there- 
fore turn our attention mainly to the quadratic term (p 2 ) 
dominant stage. Although some analytic solutions are 
already given in Ref. pl| , we first summarize those and 
then show numerical analysis. We then discuss a nat- 
ural quintessence scenario by this model and give some 
constraints for unknown parameters such as 777,5 ■ 



A. Analytic solutions in p 2 -dominant stage. 

In the conventional universe, matter energy density is 
negligibly small in comparison with radiation energy den- 
sity in the early stage of the universe. Then we know that 
it is also the case in the ^-dominant stage. Hence we con- 
sider only radiation fluid and a scalar field 4>- As for the 
initial conditions, we have two possible cases: one is the 
radiation dominant case and the other is the 0-dominant 
one. We discuss these separately. 

(1) radiation dominant initial conditions 

In this case, the scale factor e xpa nds as a a t x ' A and 
the equation for the scalar field (2.4) is now 



+ -4>-ap^ ( j ) 
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We find an analytic solution for a < 6 that is 

2 

t ^ °+ 2 



(3.3) 
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with = — 7. A* to ( 3 - 4 ) 
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wher e t$ and 4>q are integration constants. The constraint 
( |3.4|) requires a < 6. The energy density of the scalar 
field evolves as 

3a(a + 2) / t 

P<t> = Vo — 

o — a \to 
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(3.5) 



where Vq = V(<po) and ag ■■ 
parameter f2^,, we find that 



a(t ). As for the density 



Sid, = 
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= n 
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c + 2 



(3.6) 



Vo/p T (to) because p r cx a 



where f2 

If a > 2, just contrary to the tracking solution, the 
scalar field energy decreases faster tha n th at of the radi- 
ation. This result is confirmed by Eq.( |2.9[ ) with the fact 
that the kinetic energy of a scalar field pi = 4> 2 /2 turns 
out to be larger than 4U as pJ; K) /(W) = 2a/ (6 - a). 
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If a — 2, the scalar field energy drops at the same rate 
as that of the radiation until the conventional universe 
is recovered. This is the so-called "scaling" solution. If 
a < 2, the scalar field energy decreases slower than the 
radiation energy, and eventually the scalar field domi- 
nates the radiation. 

What happens for a > 6 ? Since there is no asymp- 
totic solution for which a kinetic term balances with a 
potential term, we expect that a kinetic term dominant 
solution is obtained asymptotically. It is easy to show 
that a potential dominant (or slow rolling) condition is 
not preserved asymptotically. 

Assuming that a kinetic term is dominant, Eq. ( |3.2| ) 
is now 



0+^0 = 0, 
finding an asymptotic solution 
(f> oc i 3 . 

With this solution, wc find that 



(K) 



1 



-d) z oc t 



-3/2 



2 
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(3.7) 

(3.8) 

(3.9) 
(3.10) 



where pj" and 



denote the kinetic and the po- 
tential terms of the energy density of a scalar field, re- 
spectively. Hence, if a > 6, kinetic term dominance is 
guaranteed. This is also confirmed by numerical calcula- 
tion. 

(2) (/>-dominant initial conditions 

If a scalar field initially dominates radiation, what we 
will find in the dynamics of the scalar field? In the con- 
ventional universe, once a quintessence field dominates, 
radiation dominance will not be obtained. In the present 
model, however, a radiation dominant era is recovered as 
we will see. 

Assuming the scal ar field dominance, we find the Fried- 
mann equation (2.2) as 
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6m 5 3 
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while the equation for the scalar field (2.4) is 



'3H<fi — aji 



a+4 ±— a—1 _ 
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(3.11) 



(3.12) 



Inserting Eq. (3.11) into Eq. (3.12), we find a second or- 
der differential equation for 4>. The asymptotic behavior 
of the solution can be classified into three cases: (a) slow 
rolling (a potential term dominant) solution (a < 2), (b) 
a solution in which the potential term balances with the 
kinetic term (a — 2), and (c) a kinetic term dominant 
solution (a > 2). 



First, assuming a slow rolling condition (4> 2 <C V and 
i| < H\<j>\, \V'\), Eqs. j3Tl] ) and ( gig ) are 
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From Eqs. fl3.13| ) and fl3.14| ), we obtain 
finding an analytic solution 
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(3.14) 



2am 5 3 , 
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where to is an integration constant. With this solution, 
we have (j) 2 oc t~ x and V oc t~ a l 2 . In order to guarantee 
the slow-rolling conditions, we have to require a < 2. In 
this case, the Universe expands as 



a = flo exp{[i/ (i-<o)] (2 ~ Q)/2 } 
where a constant Hq is given by 
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Note that this solution describes an inflationary evolu- 
tion, whose expansion is weaker than the conventional 
exponential inflation, but stronger than the power-law 
type. Although this solution does not play into anything 
with a quintessence scenario, it may be interesting to 
discuss a spectrum of density perturbations for such an 
inflationary scenario. The results will be published else- 
where. 

Sec ondly, we a ssume a kinetic term dominance. Eqs. 
( |3.11 ) and (3.12) are now 
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Combining Eqs. ( glf ) and ( gig ), we have 



-(1/4to 5 3 )</> 3 , finding a solution as 
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where £o an d 4>o are integration constants. For the 
solution with + sign, the potential term decreases as 
V oc i~ Q / 2 , while the kinetic term drops as <f> 2 oc 
Hence, if a > 2, the kin etic term dominance condition is 
preserved. From Eq. ( |3.18 ), we find that the Universe 
expands as 



a (tAo) 1/6 , 



(3.21) 



and then the energy density of the scalar field decreases 
as p$ oc t^ 1 oc a~ 6 , while p r oc a~ 4 . Therefore, with this 
solution, the scalar field energy decreases faster than that 
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of radiation, and eventually radiation dominance will be 
reached. 

For the solution ( |3.20 ) with — sign, the scalar field 
evolves into as t approaches a critical value i cr large, 
climbing the potential as V — + oo. However, before reach- 
ing this critical point, the potential becomes dominant, 
and then the assumption of a kinetic dominance is no 
longer valid in this limit. With the previous analysis 
for the potential dominance, we expect that this solution 
will also reach to radiation dominant stage. This is con- 
firmed by numerical calculations. As a conclusion, the 
asymptotic behavior for the case with a > 2 is described 
by kinetic term dominance of the scalar field, followed 
eventually by radiation dominance. 

For the remaining case of a — 2, we find an analytic 
solution ]2l]] , which is a power law expansion of the Uni- 
verse 



a = a {t/t ) p , 



with 



and 



P= 6 



1 + lyn, 



= 2V2m^V/ 2 . 
The scalar field energy density evolves as 
Pj, oc t~ x cx a -1 / p . 



(3.22) 



(3.23) 



(3.24) 



(3.25) 



If fi > 40 1/6 m 5 w 1.85m 5 , p > 1, that is, we find 
a power-law inflationary solution. The power-law infla- 
tion is, however, nothing to do with quintessence, al- 
though it may be interesting in the early universe. If 
fi < 4 1 / 6 m5 I.267715, p < 1/4, and then the scalar field 
energy decreases faster than that of radiation. 
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TABLE I. The fate of a scalar field for each initial condition, 
where S and R denote scalar field dominance (p a S> pr) and radi- 
ation dominance (p r 3> pd,), respectively. The "scaling" means 
the scaling solution (CIj, =constant). The asymptotic behavior 
constant is described by "const", while the ki- 



netic dominance of pi^' 3> pT is denoted by "kinetic". "PL" is 
power-law. "DE" is decelerating power-law expansion. 



Hence radiation dominance will eventually be obtained. 
For 4 1 / 6 77i5 < [i < 40 1 / 6 m5, we can easily show that this 
solution is a global attractor and the ratio of radiation 
energy to that of the scalar field remains constant |24| . 

We summarize the obtained analytic solutions and its 
fate in Table 1. 



B. Numerical analysis in /^-dominant stage 

Now we study the dynamical property of the scalar 
field numerically and show that the above analytic solu- 
tions are really attractors. The systematic analysis of the 
dynamical properties of the scalar field in the quadratic 
dominant stage will be given elsewhere [^4| . Here we con- 
sider only the case of 2 < a < 6 because we are interested 
in quintessence. 
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FIG. 1. (Top) The time evolution of p r and P0 start- 
ing from a radiation dominant initial condition. We set 
m 5 = 2.15 x 10"" 



tial conditions, 
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we set p r 

29™ 4 
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and 



' = 1.0 X IQ-^ml, 1.0 x IQ-^ml 



and 1.0 x 10- 41 m 4 4 



for (1), (2), and (3), respectively. If p,£ is initially greater 
than the value of the attractor solution (case (1)), the kinetic 
term soon dominates and eventually the attractor solution is 
reached. If p$ is initially less than the attractor's value (case 
(3)), the potential term dominates until the solution reaches 
the attractor. Thus, for a wide range of initial conditions ((1) 
~ (3)), P4, reaches that of the attractor solution. (Bottom) 
Corresponding density parameter of the scalar field (f2</,) de- 
creases after the attractor solution is reached. 
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FIG. 2. (Top) The time evolution of p r and p^ from the 
scalar-field dominant initial condition. The kinetic term soon 
dominates the energy density of the scalar field, and radia- 
tion dominance is eventually reached. As for initial condi- 
tions, we set p r = 1.0 x lO -25 ™/, p£ K) = 9.35 x 10" 20 m 4 4 , 

p^ P) = 1.0 x 10~ 17 ra 4 4 , (Bottom) The time evolution of the 
ratio of the potential term of the energy density of a scalar 
field to the total total energy density of the scalar field. This 
suggests that the kinetic term dominates the potential term 
soon. We set ms = 2.15 x 10 _3 m4, and p = 1.0 x 10 _s m4. 



First, we show that the solution (3.3) is a unique at- 
tractor for plausible initial conditions. We depict the 
evolution of p r and p^ for various initial conditions in 
Figs|] and | for a = 3. The fi gures in Fig.|l] show that 
for a wide range of initial conditions, the energy density 
of the scalar field approaches that of the attractor solu- 
tion (3.3). We also find that even starting from scalar 



field dominance, the universe eventually evolves into the 
radiation dominant stage (see Fig. 2). In fact, the ki- 
netic term of a scalar field always turns out to be domi- 
nant even if we start with a potential dominance. Since 
p^ oc a~ 6 in the case of the kinetic term dominance and 
p r oc a -4 , radiation energy eventually overcomes that of 
the scalar field and the Universe evolves into a radiation 
dominant era. As a result, the attractor solution is al- 
ways reached for any initial conditions. We find similar 
results for any value of a in 2 < a < 6. We may con- 
clude that the solution (3.3) is a unique attractor in the 
present dynamical system. 

Once the attractor is reached, the scalar field energy 
decreases faster than that of radiation, which is a most in- 
teresting feature in the brane quintessence scenario. It is 
worthwhile noting that this potential in the conventional 
cosmology without the quadratic term, the quintessence 
scenario does not work if the scalar field energy initially 
dominates that of radiation ||. Next, we shall discuss a 
more natural quintessence scenario in the brane world. 



C. Quintessence scenario 

Now we are ready to discuss a quintessence scenario in 
the brane world pifl . We assume 2 < a < 6. Using two 
attractor solutions (one in the p 2 -dominant stage and the 
other in the conventional universe), we show a successful 
and natural scenario. Since the quintessence solution in 
the conventional universe model is an attractor, our so- 
lution should also recover the same trajectory after the 
quadratic term decreases to be very small. We have con- 
firmed this numerically. The main difference is that we 
can include not only radiation dominant initial condi- 
tions but also scalar-field dominant initial conditions for 
a successful scenario. In the numerical analysis, to evalu- 
ate the present value of the density parameter of a scalar 
field, we include the matter fluid as well as radiation and 
scalar field. 

(1) a scenario 

First, we shall overview a quintessence scenario us- 
ing attractor solutions . We introduce t s (the cosmic 
time when the attractor solution in p 2 dominant stage is 
reached), t c (when the p 2 -term drops just below the con- 
ventional density term), iNs (nucleosynthesis), i cq (whcn 
radiation energy density becomes equal to matter den- 
sity), idoc(the decoupling time) and to (the present time). 
If we approximate the evolution of the Universe by the 
attractor solutions in each stage, we find the analytic 
solution for the scalar field as follows. Normalizing the 
variables by 4-dimensional Planck mass scale 777,4, we find 
that the energy density p^ in each stage is described very 
simply as 



P4> r i 



2f 



F(a) 



- 

\1ri4 



(m 4 i)~ 



(3.26) 



G 



where F(a) is a dimensionless constant denned only by 



In the /^-dominant stage, 



F(a) = {a + 2) 



6 



(3.27) 



and in the conventional universe, 
F(a) = (5a 



3c + 2 

2^+5- (a 



12) [2(a + 6)P~ 
(radiation dominant era) 
2)(a + 4)"^T2 
(matter dominant era). 



(3.28) 
(3.29) 



Since the attractor solutions are independent, when the 
Universe shifts from one attractor solution to the other 
one, we expect a discrepancy in the energy density. How- 
ever, since the difference in p^ in each stage appears only 
in the factor F(a), the discrepancy between two attrac- 
tor solutions is given by a. We can easily check that 
the ratio of F at the p 2 -dominance to that in the con- 
ventional radiation dominance is about 0.5 to 1 unless 
a w 6. Note that the ratio of radiation dominant case 
to matter dominant case is about 0.8 for any values of 
a. Hence, when t = t c , there is a little discrepancy be- 
tween the scalar field energy densities estimated by two 
attractor solutions. 

Although the energy density of a scalar field changes 
quite similarly in any stages, the radiation energy shows a 
big difference between p 2 -dominant stage and the conven- 
tional universe. In fact, the radiation density decreases 
as a -4 , but the scale factor changes as a cx t 1 / 4 in the p 2 - 
dominant stage in contrast with a cx t 1 / 2 in the conven- 
tional radiation dominant era. Therefore, when we dis- 
cuss the density parameter of a scalar field fl^ ~ p^/pr, 
its behavior in the p 2 -dominant stage is completely dif- 
ferent from that in the conventional universe. 

^ cx i( 2 - Q )/( 2 +«) ~ a ( 2 -«)/( 2 +«) (p 2 -dominant) 



cx i 4 /( 2+a > 
cx i 4 /( 2+a > 



a 



ft 



8/(2+a) 
6/(2+a) 



(radiation dominant) 
(matter dominant) (3.30) 



Since the radiation energy must be continuous, if we ig- 
nore the above small discrepancies at t c and i cq , we can 
estimate the density parameter as 



= n} s) x 



a— 2 
" 2a + 2 
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with the above analytic attractor solutions, we cannot 
properly treat the transition between p 2 -dominant stage 
and the conventional universe. We show one numerical 
result for a = 5. We set p = 6.0 x 10 _14 m4 , 7715 
fi.O x 10 _14 m4for a successful quintessence. 

As initial conditions at a., = 1, we have chosen the at- 

,2 -do min ant st age, and 



tractor solution Eq.(3.5) i n th e p" 

solve the basic equations ([2.2]), (2. 4), (2.6) and (2/7) in- 



cluding radiation and matter fluids. In Fig.^, we depict 
the energy densities of a scalar field, radiation and matter 
in terms of a scale factor a. It turns out, as we expected, 
the discrepancy at a c is very small and the evolution 
approximately follows the attractors in each stage (at- 
tractor solutions as references). 
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= 6.0 x 10 



ni4. 



(3 



( s) 

where is the density parameter when the attractor 
solution is reached. For a successful quintessence sce- 
nario, we require that the present value of the density 
parameter of the scalar field is ~ 0.7. 

Before finding a constraint, it may be useful to confirm 
the above analysis by numerical study. This is because 



6 6.5 7 7.5 8 8.5 9 

log a 

FIG. 3. (Top) The evolution of the energy densities in 
terms of a scale factor. We set p = 6.0 x 10 _14 m 4 , and 
As for the initial conditions, we 
/ P) = 2.50 x 10 - 56 m 4 4 , 
and p T = 3.78 x 10 _S4 m 4 4 . The matter fluid is also in- 
cluded to find the present matter dominant universe as 
p m = 5.43 x 10~ 69 m 4 4 . The amount of matter fluid is chosen 
in order to find the present universe. 
31)JBottom) The enlargement of the top figure around the tran- 
sition era from p 2 -dominant stage to the conventional uni- 
verse. The thin solid lines are attractor solutions both in 
p 2 -dominant stage and in the conventional universe. 



= 3.82 x 10 - 55 m 4 4 , 
3.78 x 10 _54 m 4 4 . 



7 



-e- 



10 " 7 



15 a eq 



ao 



log a 



FIG. 4. Time evolution of the density parameter of a scalar 
field for the same model in Fig.^. 



From FigM we find that 

10 -H2 



10 



-90„ 



at decoupling time, and ~ 10 



-120 



at tc . 

m£ at 



present (see Fig.|3j). These small values guarantee a suc- 
cessful quintessence. The conventional quintessence sce- 
nario (a tracking solution) is really recovered when the 
quadratic energy density becomes small enough. We also 
show the evolution of the density parameter of the scalar 
field in Fig]] 

(2) constraints to the extra-dimension 

We discuss constraints for a natural and successful 
quintessence. We consider three constraints: nucleosyn- 
thesis, matter dominance at decoupling time, and natural 
initial conditions, in order. 

(a) nucleosynthesis 

One of the most successful results of the Big-Bang stan- 
dard cosmology is a natural explanation of the present 
amount of light elements. Therefore, in any cosmologi- 
cal models, a successful nucleosynthesis provides a nec- 
essary constraint, which may be most stringent. During 
nucleosynthesis, the universe must expand as the con- 
ventional radiation dominant era. Therefore the tran- 
sition from /9 2 -dominant stage to the conventional uni- 
verse must take place before nucleosynthesis. This con- 
straint gives a lower bound for the value of ms- Intro- 
ducing two temperatures, T c and Tns, which correspond 
to those at the transition time t c and at nucleosynthesis 
time £ns , respectively, we describe the present constraint 
as T c > TnSj which implies 



Pc > Pr{t c ) = ^Sns^n!, 



(3.32) 



where <?ns is the degree of freedom of particles at nucle- 
osynthesis. Sinc e p c — \2m\/m\ and Tns ~ lMeV, the 
constraint (§]32|) yields 

m 5 > 1.6 x 10 4 (.g NS /100) 1/6 

x (T NS /1 MeV) 2/3 GeV. (3.33) 

In the Planck unit, this constraint can be written as 777,5 > 



10- 



7774. 



If the Randall-Sundrum II model is a funda- 



mental theory, in order to recover the Newtonian force 
above 1mm scale in the brane word, the 5-dimensional 
Planck mass is constrained as 7775 > 10 8 GeV ~ 4 x 
10~ n 7774 fyfl , which would be a stronger constraint. 
However, the Randall-Sundrum II model could be an ef- 
fective theory, derived from more fundamental higher- 
dimensional theories such as Hofava and E. Witten the- 
ory ]13|] . Thus, we adopt the above constraint here. 

(b) The decoupling 

From the observation of the cosmic microwave back- 
ground (CMB), we have information of the Universe at 
T ~ 4000K, from which we expect that inhomogeneity 
of the Universe was about 10~ 5 . In order to form some 
structure from the decoupling time to the present, the en- 
ergy density of matter fluid should be larger than "dark 
energy" (that of a scalar field) by a few orders of magni- 
tude at the decoupling time. 

In p 2 -dominant stage, the Friedmann equation in the 
radiation dominant era is given by 



H = 



1 



6777,5 



Pr, 



and then using a oc i 1 / 4 , we find 
3 



(3.34) 



(3.35) 



From, p c = 12777, 5 6 /m 4 2 , we find t c = m 4 2 /(8777, 5 3 ). In the 
conventional universe, p t oc t~ 2 . Then we have 



3 

16 



(3.36) 



The energy density of matter fluid at decoupling time 
is given as 



Pm(^dcc) — Pm(^eq) 



9eqT eq X 



T, 



PAtcq) 
3 



Ta, 



CC| 



(3.37) 



As for the ene rgy d ensity of a s calar field, assuming that 
the attractor (3.26) with (3.29) is reached and using Eq. 
( |3.36 ), we can estimate p^(idcc) in terms of p,T cq and 
Tdcc for each a. 

T cq is about 10 4 K < T cq < 10 5 K. In order to im- 
pose the most stringent constraint, we adopt T en = 10 4 



K here. Setting T doc = 4000 K and T t 



cq 



„ 6q 

10 4 K, from 

the constraint of p m (tdcc) > P4>(tdec), we find the upper 
bound for the value of p, i.e. 



fi < 1.23 x 10~ ib TO 4 for a 
p < 1.31 x 10~ 14 777 4 for a 
p < 2.40 x 10 _13 777 4 for a = 5. 



3, 
4. 



(3.38) 
(3.39) 
(3.40) 



The value of p is fixed if a scalar field dominates now 
(Cl^(to) ~ 0.7). However, since we do not know the value 



8 



of (i from the viewpoint of particle physics, we shall let 
its value be free. We may discuss the naturalness of the 
present model, if we do not see the coincidence problem. 

(c) Initial condition 

About initial conditions, since a quintessence solution 
is an attractor, we may not need to worry. In fact, the 
conventional quintessence will be recovered even in the 
present model. What may be better in the present model 
is that a basin of the attractor becomes larger. In par- 
ticular, the conventional quintessence will not work if a 
scalar field dominates initially, but it will still work in 
the present model. 

Nevertheless, here we will study about natural initial 
conditions in the present brane scenario. Since we assume 
that a quintessence field <f> is confined on the brane, all 
energy scales including its potential should be smaller 
than the 5-dimensional Planck scale m 5 ; 



(3.41) 



The maximally possible energy density of radiation is also 
about m 5 4 . 

With the constraints (a) and (b), we restrict two un- 
known parameters; ms and /i. In Fig. (|^), we depict 
these constraints by three solid lines in the -/i pa- 
rameter space for a — 4 and 5. 
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FIG. 5. Constraints in the /j,-m 5 parameter space. Three 
solid lines are from (a) nucleosynthesis, (b) matter dominance 
at the decoupling time, and (c) the energy scale smaller than 
the 5-D Planck mass. The observation (f2^ ~ 0.7) fixes the 
value of fi (/iobs), which is given by the vertical dot-dashed 
lines. If we assume that the potential energy of a scalar field 
is initially not very small compared with its kinetic energy, 
a successful quintessence is possible for a narrow range of 
parameters, which is shown by the shaded region. 

Then if [i is fixed to find a scalar field dominance 
right now, the 5-dimcnsional Planck mass m% is lim- 
ited. For example, for a = 5, if fx < 10~ 14 m4, the con- 
straint for ins is just coming from nucleosynthesis, while 
if [A > 10~ 14 m4, m5 > /i, which is stronger than that of 
nucleosynthesis . 

Here we invoke a further constraint which could be de- 
rived from natural initial conditions. What would be the 
natural initial conditions for a scalar field? One plausible 
condition is an equipartition of each energy density. In 
this case, the radiation energy is larger than that of the 
scalar field because the degree of freedom of all particles 
g is larger than that of the scalar field. How about the 
ratio of the kinetic energy to total one of the scalar field ? 
In the conventional quintessence scenario, the potential 
energy should be initially much smaller than the kinetic 
one. In the present model, it is not the case. This is 
because the attractor solution in the p -dominant stage 
reduces the density parameter of the scalar field. There- 
fore, we may impose natural initial conditions for a scalar 
field. To be more concrete, we focus on the potential 
by fermion condensation. After our 3-dimensional brane 
world is created, a fermion pair is condensed by a symme- 
try breaking mechanism and it behaves as a scalar filed 
with a potential V(<j>). In this case, we expect that the 
potential term should play an important role from the 
initial stage. 

We consider the cases for a = 3,4, and 5. The ini- 
tial conditions for a scalar field are classified into the 
following three cases. First, if the kinetic and potential 
energy of a scalar field are the same order of magnitude, 
the attractor solution is reached soon, and then we ex- 
pect Qy ~ 0(g s _1 ), where g s is a degree of freedom of 
particles at a s . Secondly, if the potential energy is dom- 
inant, it will not change so much before reaching the at- 
tractor solution, and then we expect that ^^>0((7 ;s ~ 1 ). 
Thirdly, if the kinetic energy is larger than the poten- 
tial one, it will decay soon, finding an attractor solution, 
then ^^<0((7 S ~ 1 ), unless the kinetic energy dominates 
a lot, which we do not assume here. Therefore, a natural 



fi/lTl^ initial condition predicts ~ 0(g s 1 ) 



From Eqs. (3.26) and (3.35), we find 



n 



(s) 



G(a) 



2(a+4.) 

joT5J 



{m 5 t s y 



(3.42) 



where G(a) 



a+4 

2"+ 2 a" 



i (a + 2)~ — (6 - a)~^+2/3. 
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Since Pi(t s )<mf, we have a constraint that m5i s >1.5 
from Eq. d3.35| ). This constraint with f^ s) ~ 0(ff a _1 ) 
and ~ 10 3 gives the lower bound for /i/7715 as 



/x>0.131m5 for a = 3, 
/i>0.140m 5 for a = 4, 
/z>0.146m5 for a = 5. 



(3.43) 
(3.44) 
(3.45) 



With the previous constraint fi < , we find a narrow 
strip in the p-m^ parameter space, which is shown by a 
shaded region in Fig.(||) The allowed region gets smaller 
for smaller values of a. In particular, we find that no 
region is allowed for a < 3. Therefore, the present model 
prefers a rather large value of a. 

The values of /z bs, which explains the observed value 
of the dark energy now, are 



/^obs r 


- 6.25 x 10" 


18 m 4 


for 


ex. 


= 3, 


(3.46) 


/^obs r 


- 8.75 x 10" 


16 m 4 


for 


a 


= 4, 


(3.47) 


/^obs r 


- 4.06 x 10" 


14 m 4 


for 


a 


= 5, 


(3.48) 



With these values, we find that there are no natural 
ranges for a < 4. For a>5, the 5-dimensional Planck 
scale is strictly constrained from the observation because 
the allowed region is very narrow. For example, for 
a — 5, we find 



4.06 x 10~ 14 m 4 <m 5 <2.78 x 10 



-13 



m 4 . 



(3.49) 



2. If A 2 < 3(wb + 1), the late time attractor is a scalar 
field dominant solution (fi^, = 1) with = — 1 + A 2 /3. 

Case 1 is the so-called scaling solution. If it is obtained 
in a radiation dominated era, a successful nucleosynthe- 
sis is possible for A 2 > 20. However, the present obser- 
vations of a scalar field dominance (f2^ ~ 0.7) cannot be 
explained by this type of solution. Case 2 is preferred in 
the context of a quintessence scenario, but a scalar field 
behaves just as a cosmological constant in its evolution. 



Then, in order to explain Q 



(0) 



0.7, an extreme fine- 



tuning in a choice of the initial value of a scalar field or 
in a mass scale /i is required just as the case of a cosmo- 
logical constant. Therefore, some modification for this 
type of potential has been done by several authors for a 
successful quintessence JhJ, [filf . 

In the present paper, we study the effects of the 
quadratic term and see whether a natural initial condi- 
tion is found. Hence, we will not analyze each modified 
potential quintessence model, but rather study the uni- 
versal properties which are found in an exponential type 
potential. In particular, we are interested in case 2 above 
and see whether a fine-tuning is loosened by the present 
scenario. 

The organization of this section is as follows: first, as in 
the previous section, we focus on the p 2 -dominant stage, 
and we present an attractor solution which is expected to 
be found as its asymptotic behavior. This is confirmed 
by numerical analysis. Then, we discuss the possibility 
to improve a quintessence scenario by the effects of the 
quadratic term. 



IV. EXPONENTIAL POTENTIAL 



Next we investigate an exponential potential model, 



i.e. 



v{<t>) 



ji A exp 



-A- 



TO 4 



(4.1) 



which is another typical potential for a quintessence , 
P). This type of potential is often found in unified the- 
ories of fundamental interactions of particles such as su- 
pergravity theory pj| . 

Within the conventional universe, this potential shows 
an interesting property, although in itself it may not pro- 
vide a successful quintessence scenario. We first recall a 
few results §, §. 

Suppose that a spatially flat FRW universe evolves 
with a scalar field 4> and a background fluid of an equa- 
tion of state pb = wbPb- There exist just two possible 
attractor solutions, which show quite different late time 
properties, depending on the values of A and wb as fol- 
lows: 



1. For A 2 > 3{w B + 1), 



the scalar field mimics a 
barotropic fluid with = p<p/p<t, — wb, and the rela- 
tion fl<p ~ p<j,/ pb = 3(wb + 1)/A 2 holds, where is the 
density parameter of the scalar field. 



A. Analytic solutions in the p 2 -dominant stage 

Since 7715 rather than m 4 is a fundamental parameter 
in the present model, it may be natural to introduce a 
new parameter A = (ms/m 4 )A, and to write the potential 
in the form of 



V{<t>) 



4 

p exp 



-A- 



(4.2) 



In the brane world scenario, the value of A is expected 
to be of order unity, unless the potential is coming from 
other physical origin. 

We discuss two initial conditions, a radiation dominant 
initial condition and a scalar-field dominant one, in that 
order. 

(1) radiation dominant initial condition 

In the case of a radiation dominant era, the scale factor 
expa nds as a cx t 1 ^ 4 and the equation for the scalar field 
.4h is now 



4r 







exp 


\m 5 1 


m 5 



= 0. 



(4.3) 
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Since the exponential potential drops much faster than 
the inverse-power potential, unless <f> <C 7715/ A, we expect 
that the kinetic energy dominant solution is asymptoti- 
cally found as the case of the inverse-power potential with 
a > 6. Hence, assuming a kinetic energy dominant con- 
dition, we analyze the equations. Ignoring the potential 
term in Eq. (B^q), we find 



0, (4.4) 
which leads to the evolution of a scalar field as 



(j) oc t 



1/4 



(4.5) 



(K) 



t~ 3 / 2 and 



(P) 



oc 



and denote a kinetic term 



From this solution 
exp(— t 1 / 4 ) , where 

of a scalar field and a potential one, respectively. This 
behavior confirms that the above kinetic-term dominant 
solution gives an asymptotic behavior of the scalar field. 

However, for the case of (/> <C m 5 /X, in particular for 
an extremely small value of A, this potential behaves al- 
most the same as a cosmological constant, and then the 
energy density of the scalar field will soon dominate radi- 
ation, although we cannot give its critical value quantita- 
tively. However, as we will see later, it will again start to 
evolve into a radiation dominant solution as long as the 
p 2 -term dominates. If the conventional universe is recov- 
ered before reaching the radiation dominant stage, then 
the radation dominance will never be obtained because 
the scalar field dominant solution is the attractor. 

(2) (/j-dominant initial condition 

If the scalar field dominates initially, the Friedmann 
equation (jO) is 



H = 



1 



6m 5 3 



1 ■ 



2 + /j exp 



-A- 



m 5 



(4.6) 



into a radiation dominant era, just as discussed in the 
previous section. 

However, as discussed before, if A is extremely small, 
the potential does not decay so fast and then after recov- 
ering the conventional universe, a scalar field will lead to 
inflation. 



B. Numerical analysis 

In order to show that the above kinetic term dominant 
solution in a radiation dominant era is a unique attractor 
for any initial condition, we present numerical results. In 
Fig|| we depict the evolution of each energy density in 
the quadratic term dominant stage. In the top figure, 
we show the results for the initial condition of radiation 
dominance. We find that even for the initial conditions 
such that the potential term of a scalar field dominates 
the kinetic one, we eventually find a kinetic dominant so- 
lution. In the bottom, we also show the case of a scalar- 
field dominant initial condition. The universe eventually 
evolves into a radiation dominant era. Therefore, as long 
as A is of order unity, the solution obtained above is found 
asymptotically both from the radiation dominant initial 
condition and from the scalar-field dominant initial con- 
dition. We may conclude that radation dominance is a 
natural condition for the p 2 -dominant stage. 
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while the equation for the scalar field (2.4) is 
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3H<P 







exp 


\m 5 1 


m 5 



= 0. 



(4.7) 



We again expect that the kinetic term dominant solu- 
tion gives an asymptotic behavior. Assuming the kinetic 
term dominant condition, from Eqs. (4.6) and (4.7), we 
find 



1 



4m 5 3 



-0, 



(4.8) 



which is the same as Eqs. Q3.18D and ( |3.19| ), because 
the potential term does play no role. We find the same 



solution for a scalar field as before (Eq.(3.20)), and then 
the same result, i.e. this gives an asymptotic solution. 
Since radiation energy decreases slower than that of a 
massless scalar field, the universe will eventually evolve 
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FIG. 6. (Top) Time evolution of each energy density for 
A = 1.0 in the radiation dominant era. The kinetic term domi- 
nant solution is realized for any initial conditions including the 
case of potential dominance. As for the initial conditions, we 
1.0 x lO" lo m 5 4 ,p (K) = O.O,p (P) = 1.39 x 10" 11 
„ ( K ) _ nn n ( p ) - Ann ^ m-22 



set p r = 

for (1), p^> = 0.0,., 

(2), p^ = 0.0, p« 
(Bottom) Time evolution 
A = 1.0, starting from the 
tial condition. As for the 



4.22 x 10 

= 6.71 x 10 
of each energy 



~rn£, 



79 4 



for 
density 

scalar-field dominant 
initial condition, we 



Pr = 1.0 x W- ia ms,p^ K) = 0.0, p^ 1 = 1.39 x 10" 
The kinetic term dominates the energy density of a scalar field 
soon, and the radiation dominance is eventually reached. 



m£, 
for 

(3). 

for 
ini- 

set 



C. A quintessence scenario 

Although a pure exponential potential may not give 
a natural quintessence model, we shall study whether 
or not a similar mechanism for the value of in the 
p 2 -dominant stage works. If it works, it may provide a 
natural initial condition for a quintessence model based 
on an exponential potential. We discuss two cases; (a) 
A - 0(1) and (b) A - 0(1) in order. 

Case (a) : Suppose A ~ 0(1) because the poten- 
tial causes the 5-dimensional origin. In this case, as 
we discussed above, the radiation dominant universe is 
an attractor and the kinetic term dominates the poten- 
tial for a scalar field. When the universe evolves into 
the conventional expansion stage, the scalar field ap- 
proaches an attractor of a scaling solution soon, because 
A = (m4/m5)A ^> 1. In fact, we find A > 10 14 from 
the constraint on m$ by nucleosynthesis. The ratio of 
the scalar field energy to radation energy, which is fixed 
by A as fi^ ~ p^/ps = 3(ws + 1)/A 2 , turns out to be 
very small in the present model. Therefore, this does 
not provide any quintessence model. In order to remedy 
it, we may need an additional potential for quintessence 
models with exponential type potentials in the conven- 
tional universe [pd| . For example, suppose that the po- 
tential is V(4>) = A exp[— \4>/m^\ + A exp[A2</>/?ri4]. If 



A 



Pc 



which 



P2 ~ 10 m4, we find that min(y) 
gives the present small cosmological constant. However, 
an introduction of such an additional potential may break 
naturalness in the present model. 

Case (b) : If we have A - 0(1), i.e. A ~ <C 1, 
the potential may behave as a cosmological constant in 
p 2 -dominant stage. We show for the case with initial con- 
ditions that p r > p^ (FigjTj). The energy density of the 
scalar field remains constant and eventually dominates 
the radiation, leading to the inflationary stage, unless 
the conventional universe is recovered. Hence this does 
not change the conventional quintessence model with ex- 
ponential potential. If the kinetic term dominates in the 
energy density of a scalar field, however, the energy den- 
sity will decrease in time. 



C -40 
M -60 

O 



— — -J2*ation 






; scalar- 


field 































10 15 20 25 

log t/rrij 1 

FIG. 7. The time evolution of the energy densities for 
extremely small A, in which case the scalar field behaves 
as the cosmological constant. Even though starting from a 
radiation dominant initial condition, the scalar field domi- 
nates the radiation, and inflation happens. In this case, we 
take A = 1.0 x 10~ 14 . As for the initial condition, we set 



Pi 



1.0 x 10 



0.0, 



(P) 



1.50 x 10 



- 20 m 5 4 . 



We then study this case in what follows. We have 
to know the value of V(4>(a c )), which is approximately 
the present value of a cosmological constant, because the 
potential term remains almost constant as a cosmolog- 
ical constant after the universe enters the conventional 
expansion stage. The evolution for a scalar field in the 



p -dominant stage is given by Eq.(|4.5|), i.e. <f> cx a cx 1/T 



V((f>(a c )) ~ V((/)(a s ))exp [-X^/rrn] 

S T S ( m 5 



V(4>(a s ))exp 



-X- 



TO4 



1/2' 



(4.9) 



In order to find V(<p(a c )) ~ p cr , the exponent in the 
r.h.s. in Eq.(fO|) should be very large. However, since 
TO5 <C 77J4, the initial value of 4> s should be large as 



s >io^ i 

m 5 



1/2 



m 5 . 



(4.10) 



This initial condition may not be natural. 

We may conclude that for exponential type potential, a 
brane world does not improve the quintessence scenario. 



V. KINETIC ALLY DRIVEN QUINTESSENCE 

There is another type of quintessence model in the con- 
ventional universe, in which the quintessential dynamics 
is driven solely by a (non-canonical) kinetic term rather 
than by a potential term |2(|, [^7) . It is called "k- 
essence" . Here we shall study it in the context of a brane 
world. 

The model Lagrangian of a scalar field is given by 



S = 



-K ((j>)X + £(0)AT 2 



(5.1) 
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where X = \ (V</>) 2 . If we introduce a new scalar field 
by 




the action ( |5.l| ) is rewritten as 

s = J ^xV^gfi®) [-x + x 2 ] . 

Among them, the model defined by 
/($) ^ 4 -«$~ Q , 



(5.2) 



(5.3) 



(5.4) 



provides a "tracking" solution ]2q] . \i is a typical mass 
scale of the system and will fix when the scalar field dom- 
inates. We then expect a similar or more interesting fea- 
ture in the ^-dominant stage. It may provide more nat- 
ural initial conditions as the inverse power-law potential 
discussed in §3. Note that a scalar field $ has not mass 
dimension but inverse-mass dimension. 

Contrary to our expectation, however, we do not find a 
solution for which the density parameter of a scalar field 
decreases. Instead, as an attractor, we have a "tracking 
solution" in which the density parameter increases slower 
than that in the conventional universe. Furthermore, we 
find a "scaling" solution as a transient attractor in the 
radiation dominant era. 



A. Analytic solutions 



(1) Model 

Fir st, w e sho rtly explain our quintessence model with 
Eqs. ( |5.3| ) and (5^). Assuming the FRW universe model, 
the "pressure", p$, and the "energy density", of a 
quintessence scalar field $ is given by 



p$ = /i 4 - Q $-" (-X + X 2 ), 



/?$ = (-X + 3X 2 ), 



(5.5) 



(5.6) 



where X = (l/2)$ 2 . Note that in order to guarantee 
the positive energy density of the scalar field , p$ , it is 
necessary that 4> 2 > (2/3). 

The field equation corresponding to fl2.4|) is 



$(1 - 3$ 2 ) + 3H(l - <j> 2 )$ - —(2 - 3$ 2 )$ 2 = 0, (5.7) 

Since this equation has a reflection symmetry ($<-»•—$), 
we discuss only the case of $ > 0. 

(2) Analytic solutions in the conventional cosmology 

First we show two analytic solutions in the conven- 
tional cosmology. We assume the radiation dominant 



era, even though a similar solution is found in the mat- 
ter dominant era. 

Then the scale factor expands as a cx t 1 / 2 and the 
equation for the scalar field is Eq. (5.7) with H = l/2t. 

We have two analytic solutions: (a) a tracking solu- 
tion, which is exact and is an attractor, and (b) a scaling 
solution, which is approximate and a tr ansient attractor. 
(a) a tracking solution: For Eq.(5.7) with H — l/2t , 



we have an exact solution for a < 2 or for a > 3 as 



which yields its energy density as 



2(2 -a) 



2(3 



1 (2-«)/2 



3(2 -a) 



r 



cx a 



(5.8) 



(5.9) 



For a < 2, the energy density of the scalar field decreases 
slower than that of radiation, i.e. the density parameter 
ft® increases as a 2 ^ 2 ~ a \ This solution is the tracking 
solution. 

For a > 3, the energy density becomes negative, al- 
though the density parameter = |p$//? r | decreases. 
This solution may not be interesting for a quintessence 
scenario because the scalar field energy never dominates. 

(b) a scaling solution: Another interesting solution is 

found in the limi t of $ 2 ^ 1. In this case, the equation 
for a scalar field (5/7) with H = l/2t is 



$ + — $ - — $ 2 = 0. 
It 4$ 



(5.10) 



It is easy to find a solution for Eq. ( 5. 10] ), which is 

$oci 2 /( 4 - a > and $oct-( 2 -")/( 4 - Q ). (5.11) 

The coefficient is an integration constant and depends on 
the initial condition. The energy density is now 



p$ cx t 2 cx a 4 , 



(5.12) 



This is nothing but a scaling solution in a radiation dom- 
inant era. 

We easily show that this solution is an attractor in the 
present system with the approximation of <I> 2 3> 1. How- 
ever, this solution ( 5.11 ) shows that the approximation 
will be eventually broken because $ decreases. Hence, 
after this approximation becomes invalid, the universe 
evolves into the tracking solution. Note that for a < 3, 
$ increases and the approximation is always valid, al- 
though the energy density is negative. 

(3) Analytic solutions in the p 2 -dominant stage 

In the p 2 -dominant stage, assuming the radiation dom- 
inant era, i.e. the evolution of the scale factor is a cx t 1 / 4 , 



the equation for a scalar field is given by Eq. (5.7) with 
if = 1/4*. 
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The same as in the conventional universe, we have two 
analytic solutions: (a) a tracking solution, which is exact 
and an attractor, and (b) a scaling solution, which is 
approximate and a transient attractor. 

(a) a tracking solution: There is an exact solution for 



< 1 or for a > 3/2, which is 



$ = 



' 3 -2a 
3(1 -a) 



(5.13) 



which yields the energy density of a scalar field as 



4(1 - a) 



3- 2a 



n (2-a)/2 



3(1 -a) 



r 



(5.14) 



The density parameter f2$ increases as a 4 ^ 1- ") for a < 1. 
This solution is again a tracking solution, although the 
rate of increase is smaller than that in the conventional 
universe. There is no solution in which 0$ decreases. 
This is a big difference between the present model and 
the model with inverse-power potential. For the case of 
a > 3/2, p$ becomes negative. Since the positivity of the 
energy density of a scalar field is not guaranteed, it might 
be allowed if the Friedmann equation is not contradicted 
as in the case of a radiation dominant era. Although it 
might be interesting because the density parameter f2$ 
decreases, it turns out that such a solution cannot evolve 
into a tracking solution in the conventional universe p6[ . 
Hence, we do not discuss this solution in what follows, 
(b) a scaling solution: We have another solution in the 

limiting case of <I> 2 3> 1. The equation for the scalar field 



$ _l $ 



0. 



At 4$ 
has an analytic solution 

Qott 3 ^ 4 ^ and $ oc t -(l-«)/<4-«) 



and 



,-1 -4 

p$ oc i cx a , 



(5.15) 



(5.16) 



(5.17) 



whi ch is a scaling solution. Since $ decreases as 



Eq.( |5.16 ), the approximation of <£> 2 ^> 1 will be even- 



tually broken. This scaling solution is a transient attrac- 
tor. Then, a tracking solution will be finally reached for 
a < 1. Since there is no attractor solution for a > 1, 
the evolution of a scalar field may depend on the initial 
conditions. Once the universe evolves into the conven- 
tional expansion stage, however, the scalar field begins 
to approach an attractor solution. 

Next, we investigate the case of the scalar-field domi- 
nant era. The Friedmann equation (|2.2|) is given by 



H = 



1 



6m 5 3 



-X + 2>X 2 



(5.18) 



aoct* 

$n 3 

M 

where a dimensionless constant $o is given by 

3 4 

4 - a 



$ = 2— 



JL 

m 5 



(5.19) 



(5.20) 



Then the energy density of the scalar field is given by 
3 



p$ = -m\t 1 oc a 4 , 



(5.21) 



which drops as the radiation energy. Hence, once this 
solution is reached, contrary to the case of inverse-power 
law potential, the radiation never dominates the scalar 
field. Therefore, a scalar field dominant initial condition 
does not provide a successful quintessence scenario. 



B. Numerical analysis 

In order to confirm the above analysis, we shall show 
our numerical results. The initial condition is chosen such 
that the inequality 4> 2 > 2/3 is satisfied, which guar- 
antees positivity of the energy density of a scalar field. 
First, we show the case with radiation dominant initial 
conditions in Fig.|| We find that a tracking solution 
Q5.13 ) or a scaling solutions (5.16) is really an attractor 
or a transient one. 



In the limit of $ 2 3> 1, we find a power law solution as 
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(Top) Time evolution of each energy density for 
26 m . „, 1 r, ^ m- 14. 



FIG 

a = 0.5. We set fi = 1.0 X 10 
As for initial conditions, we set p r = 1.0 x 10 
p* = 1.24 xl0~ 72 m 4 4 ($ = 1.15x10 
1.96 x 10" 93 m 4 4 ($ = 1.15 x 10 3 m 4 ~\X = 



~ 52 m 4 

m 4 , 



and 



6.21 x 10" 97 m 4 4 



1 ,X = 6.67xl0 11 ), 
6.67 x 10" 1 ), 



•<!> 



1.15 x lO^m^yX = 6.67 x 10" 



-In 



and 1.95 x 10" ""m^ ($ = 1.15 x 10^m 4 _i , X = 6.67 x 10" 
for (1), (2), (3) and (4), respectively. It is easy to see that the 
tracking solution is an attractor, while the scaling solution is 
a transient attractor. (Bottom) The same figure as the top 
for a = 1. We set p = 1.0 x 10~ 20 m 4 ,m 5 = 1.0 x 10 _14 m 4 . 

1.0 x 10" 52 m 4 4 , and 
- 1 X = 6.67 x HT 1 ), 



we set p r = 
= 1.15 x 10 4 m 



As for initial conditions, 

p* = 5.77 x 10 _65 m 4 4 ($ -- 

5.77 x 10 _71 m 4 4 ($ = 1.15 x 10 10 ra 4 _1 ,X = 6.67 x 10 
for (1), (2) , respectively. There is no attractor solution in 
the p 2 -dominant stage. After the conventional cosmology is 
recovered, the scaling solution evolves into a tracking solution. 

Next, we show numerical results for the scalar field 
dominant initial conditions in Fig.^. The energy den- 
sity of the radiation never dominates that of the scalar 
field, and there exists no radiation dominant era, which 
is excluded from the constraint of nucleosynthesis. 



log t/rn^ 1 

FIG. 9. (Top) Time evolution of the energy densities 
with the scalar field dominant initial conditions for a = 1. 
As for initial conditions, we set p r = 1.0 x 10 _58 m 4 4 and 
p* = 1.12xl0~ 52 m 4 4 , ($ = 1.66xl0 19 m 4 _1 , X = 2.45xl0 13 ). 
Note that both energy densities drop as oc t^ 1 in the 
p 2 -dominant era and oc t~ 2 in the conventional universe. 
(Bottom) The change of the power of scale factor, TL — Ht. 
Even in a scalar field dominant universe, the cosmic expan- 
sion is the same as that of the radiation, i.e. a oc t 1 / 4 the 
p 2 -dominant era and t 1 ^ 2 in the conventional universe. We 
set p = 1.0 x 10~ 20 m 4 ,m 5 = 1.0 x 10~ 14 m 4 . 



C. Constraint to the model 

Now, we discuss the value of the parameter fi. For a 
successful quintessence, i.e. in order for a scalar field to 
dominate the energy density right now, we have to tune 
the value of \x. Naive estimation gives 



-X + 3X£ 



P(cr)0> 



(5.22) 



where $0 and Xq are the present values of $ and X, and 
P(cr)o is the present value of the critical density. 

The present value of the scalar field is estimated by a 
tracking solution in the matter dominant era, which is 



$ = 



2(4- 



Since Hq 



8 -3a 
2/3*0, we find 



-t. 



(5.23) 
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^cT 1 = ~H 



8 - 3a 



An 



(4 -a) 



- 2a 
0(1) 



(5.24) 



8 -3a 

for < a < 8/3. These equations with Eq. ( 5.22j ) fix /i 
and the present value of the scalar field as 



M ^ 10 (43a-48)/(4-a) [GeVL 



$ ~ 10 4d [GeV -1 ]. 



(5.25) 



(5.26) 



If the mass scale /i is the same order of magnitude as the 
five-dimensional Planck mass m 5 , we have a constraint 
on the value of a from nucleosynthesis, i.e. ms>10 4 GeV, 
that is a > 1.4. 



VI. SUMMARY AND DISCUSSION 

In this paper, we have studied quitessence models in 
a brane world scenario. We have adopted the second 
Randall-Sundrum brane scenario for a concrete model, 
although a similar result would be obtained in other 
brane world models. As a consequence of a brane embed- 
ded in the extra-dimension, the quadratic term of energy 
density appears, changing the expansion law in the early 
stage of the universe. This affects the dynamics of the 
scalar field in the quadratic-term dominant stage. We 
have then investigated three candidates for a successful 
quintessence. 

As a first candidate, we have discussed an inverse- 
power-law potential model, V = /i 4+Q 0~ Q with 2 < a < 
6. We have shown the solution in which the density pa- 
rameter of a scalar field decreases as a~ 4 ( Q_2 )/( Q+2 ) in 
the /5 2 -dominant stage. This feature provides us wider 
initial conditions for a successful quintessence. In fact, 
even if the universe is initially in a scalar-field dominant, 
it eventually evolves into the radiation dominant era in 
the p 2 -dominant stage, which guarantees a successful nu- 
cleosynthesis in the conventional universe stage. 

Although initial conditions could be arbitrary because 
the present solution is an attractor, we may have a nat- 
ural initial condition for some specific origin of a poten- 
tial such as a fermion condensation. If this is the case, 
equipartition of each energy density is more likely. As- 
suming such an equipartition, we have shown constraints 
in fj, — TO5 plane for a natural and successful quintessence 
scenario. The allowed region gets wider as a is larger, 
because for larger a, the density parameter becomes 
smaller when the conventional cosmology is recovered. 
We conclude that in order to explain naturally the ob- 
servational value of the dark energy by the present sce- 
nario, a > 4 is required. This constraint also restricts the 
value of the 5-dimensional Planck mass, e.g. for a = 5, 
4.06 x Kr 14 TO 4 <m 5 <2.78 x 10~ 13 ra 4 . 

We have also discussed an exponential potential model 
V = \i A exp(— Xc/t/rrii), although by itself it may not 



provide a successful quintessence scenario. In the five- 
dimensional brane scenario, A = Xm^/m^ is expected to 
be order unity. If that is the case, the kinetic term of 
the scalar field becomes eventually dominant for any ini- 
tial conditions, and the density parameter of the scalar 
field decreases in the quadratic-term dominant stage. In 
this case, however, A becomes too large to explain the 
present scalar field dominance. We may need unnatu- 
ral modification in the potential. On the other hand, if 
A ~ 0(1) to find a successful quintessence scenario in the 
conventional universe stage, A becomes very large. This 
provides an extremely flat potential, which behaves just 
as a cosmological constant in the p 2 -dominant stage, re- 
sulting in an inflationary expansion before reaching the 
conventional universe. After that, the radiation never 
dominates, which contradicts nucleosynthesis. Therefore, 
in both cases, we may not find a natural and successful 
quintessence scenario. 

As a third model, we have investigated a kinetic-term 
quintessence (the so-called fc-essence) model. We have 
adopted a model with an inverse-power-law potential 
in coefficients of kinetic terms. This provides a track- 
ing solution just the same as the case with an inverse- 
power-law potential. We then expect to obtain a natural 
quitessence scenario. Contrary to our expectation, how- 
ever, we do not find any solution in p 2 -dominant stage 
by which the density parameter fi$ decreases. Instead, 
we find a tracking solution in which density parameter 
increases more slowly than that in the conventional uni- 
verse. We also find a scaling solution which is a transient 
attractor. Then, if the universe starts with radiation 
dominance, the density parameter keeps constant in the 
early stage and then the universe moves to a tracking 
solution, finding a usual fc-essence in the conventional 
universe. We do not find so much advantage in a brane 
world. Only the density parameter increases more slowly 
in the p 2 -dominant stage, which provides a wider ini- 
tial condition for a sucsessful quintessence. Finally , we 
have shown the value of the parameter \x appeared in 
this model can be taken as the same order as the five- 
dimensional Planck mass scale if a > 1.4. 
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